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We show effects of the event- by- event fluctuation of the initial condi- 
tions (IC) in hydrodynamic description of high-energy nuclear collisions on 
some observables. Such IC produce not only fluctuations in observables 
but, due to their bumpy structure, several non-trivial effects appear. They 
enhance production of isotropically distributed high-p^ particles, making 
V2 smaller there. Also, they reduce vi in the forward and backward regions 
where the global matter density is smaller, so where such effects become 
more efficacious. They may also produce the so-called ridge effect in the 
two large-pr particle correlation. 

PACS numbers: 25.75.-q, 24.10.Nz, 24.60.-k, 25.75.Ld 

1. Introduction 

Hydrodynamics is one of the main tools for studying the collective flow 
in high-energy nuclear collisions. In this approach, it is assumed that, after 
a complex process involving microscopic collisions of nuclear constituents, 
at a certain early instant a hot and dense matter is formed, which would be 
in local thermal equilibrium. After this instant, the system would expand 
hydrodynamically, following the well known set of differential equations. 

The initial conditions (IC) for the hydrodynamic expansion are usually 
parametrized as smooth distributions of thermodynamic quantities and four- 
velocity. However, since our systems are small, important event-by- event 



Presented at IV Workshop on Particle Correlations and Femtoscopy (WPCF2008) , 
Krakow (Poland), September 11-14, 2008. 



CD 



2 



WPCF08 PRINTED ON JANUARY 19, 2009 



fluctuations are expected in the IC of the real collisions. Moreover, each set 
of IC should presents strongly inhomogeneous structure. Our purpose here 
is to discuss some of the effects caused by such fluctuating and bumpy IC 
on observables. 

2. Event-by-event fluctuating hydrodynamics 

The main tool for our study is called NeXSPheRIO. It is a junction of 
two computational codes: NeXus and SPheRIO. The NeXus code [1] is used 
to compute the IC. It is a microscopic model based on the Regge-Gribov 
theory and, once a pair of incident nuclei and their incident energy are cho- 
sen, it can produce, in the event-by-event basis, detailed space distributions 
of energy-momentum tensor, baryon-number, strangeness and charge den- 
sities, at a given initial time r = \/t 2 — z 2 ~ 1 fm. We show in FigQ] an 
example of such a fluctuating event, produced by NeXus event generator, 
for central Au + Au collision at 200 A GeV. As seen, the energy-density dis- 
tribution is highly irregular and in a transverse plane (left panel) it presents 
several high-density blobs, whereas in a longitudinal plane (right panel) it 
presents a baton-like structure. When averaged over many events, these 
bumps disappear completely giving smooth IC, as those commonly used in 
hydro calculations. However, this bumpy structure gives important effects 
as we will show below. The SPheRIO code [2, [3], based on Smoothed Par- 
ticle Hydrodynamics (SPH) algorithm [4j, is well suited to computing the 
evolution of such systems, so complex as the one shown in FigQj. 



tau=1fm, b=0fm, eta=0 tau=1fm, b=0fm, x=-2 




-8 -6 -4 -2 2 4 6 8 -8 -6 -4 -2 2 4 6 8 

x [fm] y [fm] 

Fig. 1. Fluctuating IC, produced by NEXUS generator, for the most central Au+Au 
collisions at 200 A GeV. Left: Energy density distribution is plotted in the rj = 
plane. Right: Corresponding plot in the x = 2 fm plane. 
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3. Effects of bumpy and fluctuating IC 

In a previous paper [5], we discussed some effects of the bumpy structure 
of IC, as shown in FigUJ on spectra and V2 coefficient. The main conse- 
quence of baton structural is a violent and cylindrically isotropic expansion 
of the batons, which produces additional isotropic high-pr particles. 

3.1. Transverse-momentum spectra 

As clearly shown in Fig O this implies that high-py part of the pt 
spectrum becomes enhanced as compared with the smooth averaged IC 
case, making it more concave and closer to data. 
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Fig. 2. Charged-particle px distributions (including those from resonance decays) 
computed in two different ways. The solid line indicates result for fluctuating IC, 
whereas the dotted line the one for the averaged IC. Data points [5] are also plotted 
for comparison. 

3.2. pt dependence 0/(^2) 

As for the anisotropy of the transverse flow, we illustrate in Figj3] that 
the elliptic-flow coefficient V2 suffers reduction as we go to high-p-p region, 
due to the additional high-pj- isotropic component included now. The result 
is closer to the available data. 

3.3. 7] dependence 0/(^2) 

As for the 77 dependence of V2 , we know that the average matter density 
decreases as \n\ increases as reflected in the r] distribution of charged parti- 
cles, so the isotropic expansion effect of the batons becomes more efficacious 
there and, therefore, i>2 suffers a considerable reduction in those regions. 




1 In that paper, we have considered only the transverse structure of IC and, then called 
it granular, but all the discussions remain valid replacing granular by baton there. 
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Fig. 3. pt dependence of (i^) in the centrality window and rj interval as indicated, 
compared with data [7]. The solid line indicates result for fluctuating IC, whereas 
the dotted one that for the averaged IC. The curves are averages over PHOBOS 
centrality sub-intervals with freeze-out temperatures as indicated. 
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Fig. 4. rj dependence of (V2) for three centrality windows. The solid lines indicate 
results for fluctuating IC, whereas the dotted lines the ones for the averaged IC. 
Data points [7] are also plotted for comparison. Tj has been taken as in Figj3] 

3.4. i>2 fluctuation 

The IC fluctuation also implies a large Vi fluctuation. In preliminary 
works [8j on Au+Au collisions at 130A GeV, we showed that this actually 
happens. The results, with QGP included, have indeed been confirmed in re- 
cent experiments^ [10] . More recent computations for Au+Au at 200^4 GeV 
|llj gave similar results. In Fig. [5j we compare the latter with those data. 



3.5. Ridge effect 

Another effect, which is produced naturally by the longitudinal baton 
structure of IC, as shown in FigfTJ is the so-called ridge phenomenon which 
has been experimentally seen in high-py nearside correlations [12]. Since 
batons in fluctuating IC which are close to the surface produce longitudinally 
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Fig. 5. <Xt) 2 /(i>2) computed for Au+Au collisions at 200 A GeV, compared with data. 
Upper: <J V2 /(v2) is given as function of the impact parameter (b) and compared 
with the STAR data [S]. Lower: the same results are expressed as function of 
participant nucleon number N p and compared with the PHOBOS data, [10] . 

correlated high-p^ particles, always in the same side of the hot matter, the 
ridge structure naturally appears. This is shown in FigEl as computed with 
NeXSPheRIO by J. Takahashi et al. [13J The ridge phenomenon has been 
discussed also in connection with glasma flux tubes [Hj . 



4. Summary 

In this paper, we presented several consequences of event-by-event fluc- 
tuating IC with baton structure, computed with NeXSPheRIO code. The 
main results are: 1) The baton structure of IC produces more concave px 
spectra, as compared with smooth IC; 2) It reduces (i^) both in the high-py 
and large- 1 r; | regions; 3) It also produces the ridge structure in the nearside 
correlations; 4) Large Vi fluctuations occurs, in good agreement with data. 
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Fig. 6. Two Iiigh-pT particle correlation, computed with event-by-event fluctuating 

ic. [13] 
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